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Abstract 

We consider a stochastic billiard in a random tube which stretches 
to infinity in the direction of the first coordinate. This random tube 
is stationary and ergodic, and also it is supposed to be in some sense 
well-behaved. The stochastic billiard can be described as follows: 
when strictly inside the tube, the particle moves straight with constant 
speed. Upon hitting the boundary, it is reflected randomly, according 
to the cosine law: the density of the outgoing direction is proportional 
to the cosine of the angle between this direction and the normal vec- 
tor. We also consider the discrete-time random walk formed by the 
particle's positions at the moments of hitting the boundary. Under 
the condition of existence of the second moment of the projected jump 
length with respect to the stationary measure for the environment seen 
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from the particle, we prove the quenched invariance principles for the 
projected trajectories of the random walk and the stochastic billiard. 
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1 Introduction 

The so-called Knudsen regime in gas dynamics describes a very dilute gas 
confined between solid walls. The gas is dilute in the sense that the mean free 
path of gas molecules, i.e., the typical distance travelled between collisions 
of gas molecules, is much larger than the typical distance between consecu- 
tive collisions of the gas molecules with the walls. Hence molecules interact 
predominantly with the walls and the interaction among themselves can be 
neglected. A typical setting where this is relevant is motion of absorbed guest 
molecules in the pores of microporous solid, where both the pore diameter 
and the typical number of guest molecules inside the pores are small. 

On molecular scale the wall-molecule interaction is usually rather compli- 
cated and very difficult to handle explicitly. Hence one resorts to a stochastic 
description in which gas molecules, from now on referred to as particles, move 
ballistically between collisions with the walls, where they interact in a ran- 
dom fashion. In the Knudsen model one assumes that particles are pointlike 
and that the kinetic energy of a particle is conserved in a collision with the 
wall, but its direction of motion changes randomly. The law of this random 
reflection is taken to be the cosine law where the outgoing direction is cosine- 
distributed relative to the surface normal in the point where the particle hits 
the wall. For a motivation of this choice, sometimes also called Lambert 
reflection law, see [12]. Notice that this dynamics implies that the incoming 
direction is not relevant and is "forgotten" once a collision has happened. 
Thus this process defines a Markov chain which we call "Knudsen stochas- 
tic billiard" (KSB). The random sequence of hitting points is referred to as 
"Knudsen random walk" (KRW) [6]. 

Pores in microporous solids may have a very complicated surface. Among 
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the many possibilities an elongated tube-shaped pore surface has recently 
attracted very considerable attention [23]. A three-dimensional connected 
network of "tubes" may be regarded as constituting the entire (non-simply 
connected) interior empty space of microporous grain in which particles can 
move. In this setting parts of the surface of the individual tubes are open and 
connect to a neighboring pore so that particles can move from one to another 
pore. It is of great interest to study the large scale motion of a molecule 
along the direction in which the tube has its greatest elongation. We think 
of the direction of longest elongation of a single tube as the first coordinate 
in d- dimensional Euclidean space. Together with the locally usually very 
complicated surface of pores this leads us to introduce the notion of a random 
tube with a random surface, to be defined precisely below. 

Knudsen motion in a tube has been studied heuristically for simple reg- 
ular geometries such as a circular pipe and also numerically for self-similar 
random surfaces, see e.g. [3 [9], [21]. For the straight infinite pipe, it is not dif- 
ficult to see that in dimensions larger than two the mean-square displacement 
grows asymptotically linearly in time, i.e., diffusively, while in two dimensions 
the motion is super diffusive due to sufficiently large probability for very long 
flights between collisions. Interestingly though, rigorous work on this concep- 
tually simple problem is rare. In fact, it is not even established under which 
conditions on the pore surface the motion of a Knudsen particle has diffusive 
mean square displacement and converges to Brownian motion. Indeed, it is 
probable that one may construct counterexamples to Brownian motion in 
three and more dimensions without having to invent physically pathologi- 
cal pore surfaces, but considering a nonstationary (expanding or shrinking) 
tube instead. We refer here to the work [i20j (there, only two-dimensional 
case is considered, but it seems reasonable that one may expect similar phe- 
nomena in the general case as well). Even for the stationary tube, it seems 
reasonable that the presence of bottlenecks with random (not bounded away 
from 0) width may cause the process to be subdiffusive. 

In this work we shall define a class of single infinite random tubes for 
which we prove convergence of KSB and KRW to Brownian motion. Together 
with related earlier and on-going work this lays the foundation for addressing 
more subtle issues that arise in the study of motion in open domains which 
are finite and which allow for injection and extraction of particles. This latter 
problem is of great importance for studying the relation between the so-called 
self-diffusion coefficient, given by the asymptotic mean square displacement 
of a particle in an infinite tube under equilibrium conditions, on the one hand. 
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and the transport diffusion coefficient on the other hand. The transport 
diffusion coefficient is given by the non-equihbrium ffux in a finite open tube 
with a density gradient between two open ends of the tube. Often only one of 
these quantities can be measured experimentally, but both may be required 
for the correct interpretation of other experimental data. Hence one would 
like to investigate whether both diffusion coefficients are equal and specify 
conditions under which this the case. This is the subject of a forthcoming 
paper [7]. Here we focus on the question of diffusion in the infinite tube in 
the absence of fluxes. 

For a description of our strategy we first come to the modeling of the infi- 
nite tube. The tube stretches to infinity in the first direction. The collection 
of its sections in the transverse direction can be thought as a "well-behaved" 
function u) : a ^ Ua of the first coordinate a which values are subsets of 
]R<^-i. We assume that the boundary of the tube is Lipschitz and almost 
everywhere different iable, in order to define the process. We also assume the 
transverse sections are bounded and that there exist no long dead ends or 
too thin bottlenecks. For our long time asymptotics of the walk, we need 
some large scale homogeneity assumptions on the tube: it is quite natural to 
assume that the process ui = {ua, a G M) is random, stationary and ergodic. 
Now, the process essentially looks like a one-dimensional random walk in a 
random environment, defined by random conductances since Knudsen ran- 
dom walk is reversible. The tube serves as a random environment for our 
walk. The tube as seen from the walker is a Markov process which has a (re- 
versible) invariant law. From this picture, we understand that the random 
medium is homogenized on large scales, and that, for almost every envi- 
ronment, the walk is asymptotically Gaussian in the longitudinal direction. 
More precisely, we will prove that after the usual rescaling the trajectory of 
the KRW converges weakly to Brownian motion with some diffusion constant 
(and from this we deduce also the analogous result for the KSB). This will be 
done by showing ergodicity for the environment seen from the particle and 
using the classical "corrector approach" adapted from [T7] . 

The point of view of the particle has revealed useful [Ul [IDl E] to study 
reversible random walks in a random environment, and to obtain the central 
limit theorem for the annealed law (i.e., in the mean with respect to the 
environment). The authors from [H] obtain an (annealed) invariance prin- 
ciple for a random walk on a random point process in the Euclidean space, 
yielding an upper bound on the effective diffusivity, which agrees with the 
predictions of Mott's law. The corrector approach, by correcting the walk 
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into a martingale in a fixed environment, has been widely used to obtain 
the quenched invariance principle for the simple random walk on the infinite 
cluster in a supercritical percolation [22] in dimension (i > 4, [3] for d > 2. 
These last two references apply to walks defined by random conductances 
which are bounded and bounded away from 0. The authors in |1] and [19] 
gave a proof under the only condition of bounded conductances. 

In this paper we will leave untouched the questions of deriving heat kernel 
estimates or spectral estimates, see [T| and [13], respectively, for the corre- 
sponding results for the simple random walk on the infinite cluster. We will 
not need such estimates to show that the corrector can be neglected in the 
limit. Instead we will benefit from the essentially one-dimensional structure 
of our problem and use the Ergodic Theorem; this last ingredient will require 
introducing reference points in Section [3741 below to recover stationarity. 

The paper is organized as follows. In Section [2] we formally define the 
random tube and construct the stochastic billiard and the random walk, and 
then formulate our results. In Section [3] the process of the environment seen 
from the particle is defined and its properties are discussed. Namely, in 
Section 13.11 we define the process of the environment seen from the discrete- 
time random walk, and then prove that this process is reversible with an 
explicit reversible measure. For the sake of completeness, we do the same for 
the continuous-time stochastic billiard in Section [372| even though the results 
of this section are not needed for the subsequent discussion. In Section 13.31 
we construct the corrector function and in Section 13.41 we show that the 
corrector behaves sublinearly along a certain stationary sequence of reference 
points and that one may control the fluctuations of the corrector outside 
this sequence. Based on the machinery developed in Section [3l we give the 
proofs of our results in Section HI In Section 14.11 we prove the quenched 
invariance principle for the discrete-time random walk, and in Section 14.21 
we discuss the question of finiteness of the averaged second moment of the 
projected jump length. In Section 14.31 we prove the quenched invariance 
principle for the continuous-time stochastic billiard, also obtaining an explicit 
relation between the corresponding diffusion constants. Finally, in Appendix 
we discuss the general case of random tubes where vertical walls are allowed. 
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2 Notations and results 



Let us formally define the random tube in W^, d > 2. In this paper, R'^"^ 
will always stand for the linear subspace of which is perpendicular to the 
first coordinate vector e, we use the notation || • || for the Euclidean norm 
in M'^ or M."^'^. For k e {d-l,d} let B{x,e) = {y G : ||x - y\\ < s} be 
the open ^-neighborhood of a; G M^. Define W^^^ = G : = 1} to be 
the unit sphere in W^, and let S''"^ = E>'^~^ flM'^"^ be the unit sphere in M'^"-'-. 
We write \A\ for the /c-dimensional Lebesgue measure in case A C M*^, and 
/c-dimensional Hausdorff measure in case ^4 C S'^. Let 

Sh = {w e S>'^-^ ■.h-w>Q} 

be the half-sphere looking in the direction h. For x G M.'^, sometimes it will be 
convenient to write x = {a, u), being a the first coordinate of x and u G M*^"^; 
then, a = X ■ e, and we write u = Ux, being U the projector on W^~^. Fix 
some positive constant M, and define 

A = {h G W^-^ : ll-ull < M}. (1) 

Let A be an open domain in W^~^ or R''. We denote by dA the boundary 
of A and by A — AU dA the closure of A. 

Definition 2.1 Let k E {d — l,d}, and suppose that A is an open domain 
in M^. We say that dA is {i, L)-Lipschitz, if for any x G dA there exist an 
affine isometry Ja, : R*^ ^ R*^ and a function fx : 'MJ'~^ — > M such that 

• fx satisfies Lipschitz condition with constant L, i.e., \fx{z) — fx{z')\ < 
L\\z — z'W for all z, z' ; 

• Itx^ = 0, /a;(0) = 0, and 

3x{A n B{x, e)) ^{ze B{0, e) : ^« > /,(^«, . . . , z^'^-'^)}. 

In the degenerate case k = 1 we adopt the convention that dA is {i,L)- 
Lipschitz for any positive L if points in dA have mutual distance larger than i. 

Fix M > 0, and define for n > 1 to be the set of all open domains A 
such that A C A and dA is n)-Lipschitz. We turn ^ = U„>iC„ into 
a metric space by defining the distance between A and B to be equal to 
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\{A\B)U {B\A)\, making (S a Polish space. Let Q = C^(R) be the space of 
all continuous functions M i— €, let A be the sigma- algebra generated by the 
cylinder sets with respect to the Borel sigma- algebra on <E, and let P be a 
probability measure on {Q, A). This defines a (S- valued process u = {ua, ol G 
M) with continuous trajectories. Write da for the spatial shift: daUJ- = w.+a. 
We suppose that the process u is stationary and ergodic with respect to the 
family of shifts {6a, a G M). With a slight abuse of notation, we denote also 
by 

uj = {{a,u) G M'^ : M e oOa} 

the random domain ("tube") where the billiard lives. Intuitively, Ua is the 
"slice" obtained by crossing u with the hyperplane {a} x M.'^~^. One can 
check that, under Condition L below, the domain u is an open subset of M"^, 
and we will assume that it is connected. 
We assume the following 

Condition L. There exist e, L such that du is [e, L)-Lipschitz (in the sense 
of Definition El]) P-a.s. 

Let /i^ be the {d — 2) -dimensional Hausdorff measure on dua (in the case 
c? = 2, /i^ is simply the counting measure), and denote fi'^j^ = f^a'^ /^^- Since 
it always holds that dua C A, we can regard //^ as a measure on A (with 
supp /i^ = duJa), and fi'^j^ as a measure on A^ (with supp /x^^^ = dua x duis). 
Also, we denote by z/"^ the {d — l)-dimensional Hausdorff measure on duj; 
from Condition L one obtains that z/"^ is locally finite. 

We keep the usual notation dx, dv, dh, . . . for the {d — l)-dimensional 
Lebesgue measure on A (usually restricted to Ua for some a) or the surface 
measure on E>'^~^. 

For all X = {a,u) G du where they exist, define the normal vector 
na;(x) = n^{a,u) G S'^^^ pointing inside the domain u, and the vector 
T^uj{x) = ^oj{ci,u) G S'^^^, which is the normal vector at m G dua pointing in- 
side the domain Ua (in fact, is the normalized projection of n;^ onto M'^"^), 
see Figure [H Denote also 

l^x = l^a,u = Ho; (3^) ' ^(^(x). 

Observe that k also depends on uj, but we prefer not to write it explicitly in 
order not to overload the notations. Define the set of regular points 

T^uj = {x E duj : duj is continuously differentiable in x, |ni^(x) ■ e| 7^ 1}. 
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Figure 1: Notations for the random tube. Note that the sections may be 
disconnected 

Note that G (0, 1] for all x G TZ^^. Clearly, it holds that 

dv^{a,u) = K-ld^^l{u)da, (2) 

see Figure [21 

We suppose that the following condition holds: 

Condition R. We have v^{du \ 7^^) = 0, P-a.s. 

We say that ?/ G a) is seen from x G if there exists h G 'B'^^^ and to > 
such that X + th & oj ioT all t G (0, to) and x + toh = y. Clearly, if y is seen 
from X then x is seen from y, and we write "x ^ y" when this occurs. 

One of the main objects of study in this paper is the Knudsen random 
walk (KRW) ^ = (^n)neN, which is a discrete time Markov process on du, 
cf. [6j. It is defined through its transition density K: for x,y E duj 

i^i N ((y-x) ■n^(x))((x-y) ■n^(y)) - 1 /ON 

K[x,y) = ^d- II G 7^<^,x ^ y}, (3) 

11-^ y 1 1 

where 7^ = {j h ■ edh) ^ is the normalizing constant. This means that, 

being Puj,E^ the quenched (i.e., with fixed a;) probability and expectation, 
for any x G TZ^^ and any measurable B G duj we have 

P^[e„+i eB\^„ = x] = j K{x,y) dv^{y). 

B 
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du^{x) 




a 



da 



Figure 2: On the formula ([2]); note that 



COS(p 



Following |6], we shortly explain why this Markov chain is of natural interest. 
From C,n = X, the next step = y is performed by picking randomly the 
direction h = {y — x) / \\y — x\\ of the step according to Knudsen's cosine 
density 7^ h ■ n^^ix) dh on the half unit-sphere looking towards the interior 
of the domain. By elementary geometric considerations, one can check that 
dv^ijj) = {h ■ Yii^{ii))~^\\y — xW^'^dh and recover the previous formulas. 
Let us define also 



From ([3]) we see that K{-, ■) is symmetric, that is, K{x,y) = K{y,x) for 
all x,y E 71^; consequently, $ has this property as well: 



Clearly, both K and $ depend on uj as well, but we usually do not indicate 
this in the notations in order to keep them simple. When we have to do it, 
we write K'^, instead oi K,^. For any 7 we have 



(4) 



u, P, v) = V, a, u) 



for all a, /9 G M, M G duja, v G dup. (5) 



^'^-'^{a, u, p, v) = + 7, M, + 7, t;). 



(6) 



Moreover, the symmetry implies that 



ir-((0,M),y) =K^--((0,Wy),(-ye,n)), 



(7) 
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$'^(0, u, a, v) = $^'^'^(0, V, -a, u). 



(8) 



We need also to assume the following technical condition: 

Condition P. There exist constants N, e, S such that for P-almost every u, 
for any x,y E TZ^j with \ {x — y)-e\ < 2 there exist Bi, . . . ,Bn C du, n < N — 1 
with u'^{Bi) > 6 for alH = 1, . . . , n, and such that 

• K{x, z) > e for all z E -Bi, 

• K{y, z) > e for all z G 

• K{z, z') > e for all z G Bi, z' G -Bj+i, i = 1, . . . ,n — 1 

(if = 1 we only require that K{x,y) > e). In other words, there exists a 
"thick" path of length at most joining x and y. 

Now, following [6J, we define also the Knudsen stochastic billiard (KSB) 
{X, V). First, we do that for the process starting on the boundary dut from 
the point xq G 71^ C du. Let xq = ^0,^1,^2,^3, • • • be the trajectory of the 
random walk, and define 

n 

fc=l 

Then, for t G [T„,r„_|_i), define 

t — Tn 

IIU+1 - Ull 

The quantity Xf stands for the position of the particle at time t. Since (Xf)f>o 
is not a Markov process by itself, we define also the cadlag version of the 
motion direction at time t: 

elO e 

Then, Vt G S'^^^ and the couple (Xj, Vt)t>o is a Markov process. Of course, 
we can define also the stochastic billiard starting from the interior of uj by 
specifying its initial position Xq and initial direction Vq. 
Define 

& = {{uj.u) : uj E VL,u E duJo). 
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One of the most important objects in this paper is the probabihty measure Q 
on & defined by 

dQico, u) = ^Ko'u dpi^H dP{u;), (9) 

where Z — J^dF Jj^K,Ql^dnQ{u) is the normahzing constant. (We will show 
that Q is the invariant law of the environment seen from the walker.) To 
see that Z is finite, note that Z = f^dP da fj^'^ali^f^ai'^) translation 
invariance, that is the expected surface area of the tube restricted to the 
slab [0, 1] X M*^"^, which is finite by Condition L. Let -^^(6) be the space of 
Q-square integrable functions / : © i— > R. We use the notation (/) for the 
Q-expectation of /: 

^hJ'^^J dl4{u)K^^J{u,u), 
n A 

and we define the scalar product (•, •)^ in -^^((S) by 

(/,^>^ -^jd^j ^/^oN K-ol,f{u;,u)g{u,u). (10) 
n A 

Note that (/)^ = (l, /)^, where l(a;, u) — 1 for all u, u. 

Now, for {j3, u) G TZuj we define the local drift and the second moment of 
the jump projected on the horizontal direction: 

i^p{uj,u) = E,((6 - ^o) • e I ^0 = {(i,u)) 

{x-e- l3)K{{p,u),x) diy'^ix) 



+0O 



= J {a -P) da J diJ^^{v) Kp^u^{l3,u,a,v), (11) 

— oo A 

6^(w,«) = e;,((K, - i„) ■ e.f I e„ = (/?,«)) 



+0O 



j {a-Pfda J dji'^iv) Kp^u^{l3,u,a,v). 
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When /3 = 0, we write simply A{u,u) and b{u!,u) instead of Aq{u!,u) and 
6o(c^,w)- In Section [3] we show that (^)^ = (see (^^)- 

Let Z^^^ be the polygonal interpolation oin/m m~^/^^„ ■ e. Our main 
result is the quenched invariance principle for the horizontal projection of 
the random walk: 

Theorem 2.1 Assume Conditions L, P, R, and suppose that 

{\ < ^- (12) 

Then, there exists a constant a > such that for F- almost all uj, for any 
starting point from TZ^, a~^Z^"^^ converges in law, under P^, to Brownian 
motion as m ^ oo. 

The constant a is defined by ( 147|) below. Next, we obtain the corre- 
sponding result for the continuous time Knudsen stochastic billiard. Define 
zj:^'' = s^'^^'^Xst ■ e. Recall also a notation from [6j: for x G cj, f G S*^"^, 
define (with the convention inf = oo) 

hx{v) = X + f inf{t > : X + tv E duj} G {duj, oo}, 

so that hx{v) is the next point where the particle hits the boundary when 
starting at the location x with the direction v. 

Theorem 2.2 Assume Conditions L, P, R, and suppose that f|T2l) holds. 
Denote 



where a is from Theorem \2.1\ Then, for¥-almost all uj, for any initial condi- 
tions {xo,Vq) such that h^^lvo) G TZ^^, a'^z')^^ converges in law to Brownian 
motion as s ^ oo. 

Next, we discuss the question of validity of (fT2|) . 

Proposition 2.1 If d > 3 then ^^ holds. 



\i d = 2, then one cannot expect (U21) to be valid in the situation when u 
contains an infinite straight cylinder. Indeed, we have the following 

Proposition 2.2 In the two-dimensional case, suppose that there exists an 
interval 5 C A such that R x S* C for ¥-a.a. uj. Then (6) = oo. 
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On the other hand, with Ra{uj,u) = supda — /?|; {P,v) («,«)}, it is 
clear that (fT^ holds when Rq{(jj^u) < const for all u G ujq, P-a.s. Such an 
example is given by the tube {(a, m) G : cos a <u< cosq; + 1}, a random 
shift to make it stationary and ergodic (but not mixing). 

Remark 2.1 (i) The continuity assumption of the map a ^ uja has a geo- 
metric meaning: it prevents the tube from having "vertical walls" of nonzero 
surface measure. The reader may wonder what happens without it. First, 
the disintegration formula ([2]) of the surface measure on duj becomes a 
product dfi'^{u) d(f^{a) , where /i^ is a measure on the section of du by the 
vertical hyperplane at a, and where d(f^{a) = K,~]j^da+d(f'^{a) with a singular 
part (p'^ . If the singular part has atoms, one can see that the invariant law Q 
(see above) of the environment seen from the particle has a marginal in uo 
which is singular with respect to P. This happens because, if the vertical walls 
constitute a positive proportion of the tube 's surface, in the equilibrium the 
particle finds itself on a vertical wall with positive probability; on the other 
hand, if u) has the law P, a.s. there is no vertical wall at the origin. The gen- 
eral situation is interesting but complicated; in any case, our results continue 
to be valid in this situation as well (an important observation is that fH21) 
below would still hold, possibly with another constant). To keep things sim- 
ple, we will consider only, all through the paper, random tubes satisfying the 
continuity assumption. In the Appendix, we discuss the general case in more 
detail. Another possible approach to this general case is to work with the 
continuous-time stochastic billiard directly, cf. Section U7B . 

(a) A particular example of tubes is given by rotation invariant tubes. 
They are obtained by rotating around the first axis the graph of a positive 
bounded function. The main simplification is that, with the proper formalism, 
one can forget the transverse component u. Then the problem becomes purely 
one-dimensional. 
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3 Environment viewed from the particle and 
the construction of the corrector 



3.1 Environment viewed from the particle: discrete 
case 

One of the main methods we use in this paper is considering the environ- 
ment uj seen from the current location of the random walk. The "environment 
viewed from the particle" is the Markov chain 

((%„.eo;,Wa), ^ = 0,1,2,...) 

with state space &. The transition operator G for this process acts on func- 
tions / : e ^ R as follows (cf. (E]) and (HD) 

Gf{uj,u) = E^(/(%.eCU,W6) I eo = (0,m)) 

= j K{{0,u),x)f{e,.,u;,Ux) du^ix) 

du) 

+ 00 

= j da j djji^iv) KQ^uf{daUJ,v)^{Q,u,a,v). (13) 

-co A 

Remark 3.1 Note that our environment consists not only of the tube with an 
appropriate horizontal shift, hut also of the transverse component of the walk. 
Another possible choice for the environment would he obtained by rotating the 
shifted tube to make it pass through the origin with inner normal at this point 
given by the last coordinate vector. However, we made the present choice to 
keep notations simple. 

Next, we show that this new Markov chain is reversible with reversible 
measure Q given by (Q, which means that G is a self-adjoint operator in 
L2(6)=L2(6,Q): 

Lemma 3.1 For all f,g E -^^(S) we have (^g,Gf')^ = (^f,Gg')^. Hence, the 
law Q is invariant for the Markov chain of the environment viewed from the 
particle, which means that for any f E L^{&) and all n 

(EMiO^.-eUjMn) I eo = {0,U)]) = (/> (14) 
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Proof. Indeed, from and (fT^ . 

+00 

{9,Gf\=^ j dF j dii%{u)g{u,uy-,l^ j da 

n A -00 



1 



A 

+00 



da J d¥ J diJ,Q^^{u,v) g{uj,u)f{6aid,v)(^{0,u,a,v) (15) 

-00 n A2 
+00 

j da j dF j dfi'_%iu,v)gico,u)fie^u,v)<!>'-'^{-a,u,0,v) 
-00 n A2 

(16) 

+00 

--^ j da J dF j dfi':^^Q{u,v) g{9_a'^,u)f{LJ,v)^{-a,u,0,v) 

(17) 



-00 n A2 



1 



rfP y rf/i^(t;)/(^,t;)/s:o^i y rfa 
ri A -co 

X j dij!^^{u) Ko^y^{0,v,-a,u)g{9^a(^,u), (18) 

where we used ([6]) to pass from (|T5|) to f|T6l) . the translation invariance of P to 
pass from f|T6l) to f|T7j) . the symmetry property ([5]) to pass from f|T71) to f[T8|) . 
and the change of variable a t-^ —a to obtain the last line. □ 

Let us define a semi-definite scalar product {g,f)^ '■= {g,{,I — 
Again using (fTSjl . the translation invariance of P and the symmetry of $, we 
obtain 

+00 

{9,1}^=^ j ^'^ J ^^J dfJ'o,ai^^v)^i'^^u^Oi,v)g{uj,u){f{uj,u) - f{9aUJ,v)) 
-00 n A2 

15 



+ 00 

1 



da J d¥ J djj,'^^ Q{u,v) 

oo n A2 



X $(-«, u, 0, v)g{e^^uj, u){f{e_aUJ, u) - f{uj, f )) 



hoo 



1 

z 



da J dF J dij,'^^^{u,v) 

oo n A2 



X $(0, u, a, v)g{eaUJ, t;)(/(6'„u;, v) - f{uj, u)). 

Consequently, 

{9,f \ =^ / / / ^/"o,a(w>^)'^(0,M,a,^^) 

-oo Q A2 

X (/(^^, ^) - fidaUJ, v)){g{uj, u) - g{9aUJ, v)), 
so the Dirichlet form can be explicitly written as 

+ 00 

1 



{fj\ = ^ I da I dF I df,'^,ju,v)^0,u,a,v){f{uj,u)-f{e^uj,v)y, 

(19) 



-oo n A2 



n A 



2Z 

or, by ([2]) and 

X J diy'^{x)K{{0,u),x){f{u;,u)-f{e^.eUjMx))^. (20) 

At this point it is convenient to establish the following result: 

Lemma 3.2 T/ie Markov process with initial law Q anc? transition opera- 
tor G is ergodic. 

Proof. Suppose that / G -^^^(6) is such that / = Gf. Then <^/, = and 
so, by the translation invariance and fl20l) . 

^^dfi^iu) j du'^{x)K{{s,u),x){f{9su;,u) - f{9,.,u;,Ux)f = 

n A duj 

16 



for any s. Integrating the above equation in s and using ([2]), we obtain 
jdF j dv''{x)dv-{y)K{x,y){f{e.,.,uMx)~ f{dy.eUjMy)Y = (21) 

n (aa;)2 

Let us recall Lemma 3.3 (iii) from [6]: if for some x,y E IZ^^ we have K{x^ y) > 
0, then there exist 5 > and two neighborhoods of x and By of y such 
that K{x',y') > S for all x' G Bx,y' G By. Now, for such x,y, (pT]) implies 
that there exists a constant c{Lj,x,y) such that f{9z.e^,Uz) = c{uj,x,y) for 
z/'^-almost all 2; G B^ U By. By the irreducibility Condition P (in fact, a 
much weaker irreducibility condition would be sufficient), this implies that 
f{6z.e^Mz) = c{ijj) for z/^^-almost all z G TZ^j- Since P is ergodic, this means 
that /(cu, u) = c for /ig-almost all u and P-almost all uj. □ 



3.2 Environment viewed from the particle: continuous 
case 

For the sake of completeness, we present also an analogous result for the 
Knudsen stochastic billiard (Xt,V^). The notations and the results of this 
section will not be used in the rest of this paper. 

Let & = {{u,u) : G G ujq}, and let Q be the probability measure 
on (3 X E)'^^^ defined by 

dQ{u, u, h) = ^I{u G Uq} du dh dF{u), 
2 

where Z = \ uJo\dF is the normalizing constant. The scalar product 

(-,■)_ in L^(Q) is given, for two Q-square integrable functions f,g : 6 x 

S'^-i t-. M, by 

(/'5')Q = iy" dF j du j dhf{LJ,u,h)g{uj,u,h). 

For the continuous time KSB, the "environment viewed from the particle" is 
the Markov process {{ex,.eUJ,UXt, Vt),t> 0) with the state space 6 x S'^^K 

The transition semi-group Vt for this process acts on functions / : (5 xS"'^^ 1— >• 
M in the following way: 

VJ{uj,u,h) = E^{f{ex,.eUJ,l{Xt,Vt) I Xo = {0,u),Vo = h). 
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We show that V is quasi-reversible with respect to the law Q: 
Lemma 3.3 For all f,g& L'^iQ) t > we have 

{g,Vtf\ = {lVtg\. (22) 

with f{oi!,u,h) = /(a;,M, — /i). In particular, the law Q is invariant for the 
process {{6xfeU},KXt,Vt),t > O) . 

Proof. We first prove that fl22p implies that the law Q is invariant. Indeed, 
taking ^ = 1, we get for all test functions /, 

= ai), 

by the change of variable h into —h in the integral. Hence Q is invariant. 

We now turn to the proof of fl22l) . Introducing the notation for the 
transition kernel of KSB, 

P^{Xt G dx', Vt e dh' \ Xq = x,Vo = h) = P^{x, h; x', h')dx'dh', 

we observe that 

Vtf{u;,u,h) = J dx' J dh' P'^ {{0,u), h; x' , h') f{9^K^,uj,Ux' , h'). 

In Theorem 2.5 in [6j, it was shown that is itself quasi-reversible, i.e., 
P^(x, h; x' , h') = Pt{x', —h'; x, —h). 

Therefore, 

{9Af\ = iJdFjdu j dhg{u;,u,h)Vj{iu,u,h) 
1 

~ Z 



du j dh g{u!,u, h) J da J du j dh' 
X Pt{{0,u),h;{a,u'),h')f{e^,uj,u',h') 
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z 



— J dF J du J dh g{uj,u, h) J da J du j dh' 

X P^- ((0, u'), -h'; {-a, u), -h) u', h') 

— j da j dF j du j dhg{9_a,^,u,h) J du J dh' 

X P- ((0, u'), -h'; (-«, u), -h)f{u;, u', h') 

— j da j dF j du j dhg{9a,u},u,—h) J du J dh' 

X P^{{0,u'),h'-{a,u),h)f{uj,u',-h') 

where we used that Lebesgue measure on is product to get the second 
hne, quasi-reversibihty for the third one, Fubini and translation invariance 
of P for the fourth one, and change of variables (/i, h' , a) to (—/?-, —/?-', —a) in 
the fifth one. □ 

3.3 Construction of the corrector 

Now, we are going to construct the corrector function for the random walk i^. 
Let us show that for any g G L'^{&) 

Indeed, from ffTTl) we obtain 

+ 00 

{g,^)^ = -^ j ada j dF j fiQ^^{u,v) g{uj,u)^{0,u,a,v) 



-oo Q A2 

+ 00 



1 

z 



ada d¥ fi'^^Q{u,v) g{6_a^,u)^{—a,u,0,v) 



+ 00 



1 

z 



ada J d¥ J fiQ^^{u,v) g{9aUJ,v)^{0,u, a,v), 

-oo n A2 
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so 



A2 



(24) 

Using the Cauchy-Schwarz inequality in fl24l) . we obtain 



+00 



-00 n A2 

+00 



X— / da I dF I ij,^,^{u,v)^{0,u,a,v){g{uj,u) - g{eauj,v)) 



1/2 

2 ' 



A2 



which shows fl23|) . 

Note that, from fl2^ with g = 1 we obtain the important property 



(A) =0. 

As shown in Chapter 1 of [16], we have the variational formula 

((/ - G)-'/'A, (/ - G)-V2A)^ = (A, (/ - G)-'A)^ 

= sup{(^, A)^ - (g, g)^, (g, g)^ < 00]. 

Then, provided that f|T2l) holds, the inequality fl23|) implies that the drift A 
belongs to the range of (/ — G)^/^, and so the time- variance of A is finite. At 
this point we mention that this already implies weaker forms of the CLT, by 
applying [TJ] (under the invariant measure, or in probability with respect to 
the environment) or [TUj (under the annealed measure). With this observa- 
tion, we could have used the resolvent method originally developed in [151 [IS] 
to construct the corrector. However, it is more direct to use the method of 
the orthogonal projections on the potential subspace, cf. [H [191 [18] . 
For a; e f2, M G Ouq, define 

yt>,u = {y^duj:ye> 0,K((0,M),y) > 0}, 
= {yedij:ye<Q, K{{0, u),y) > 0}. 
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Then, in addition to the space (5, we define two spaces D'T C 971 in the 
following way: 

^= {{iu,u,y) : uj e d,u e dujo, y E V+„}, 
DJl = {{uj , u, y) : LJ E Q, u E dcuo, y G du}. 



On we define the measure KQ with mass less than 1, for which a non- 
negative function / : i-^ M has expected value 

(/>,, = ( / fiu;,u,y)K{iO,u),y)du-iy))^. (25) 

Two square-integrable functions f,g : M. have scalar product 

(^'^>K« = ( / f{^,^,y)9i^,^^,y)K{{0,u),y)du^{y)y (26) 



Also, define the gradient V as the map that transfers a function / : 6 i— >■ 
M to the function V/ : D^I ^ M, defined by 

{Vf){u,u,y) = f{9y.,u;,Uy) - fiu,u). (27) 

Since Q is reversible for G, we can write 

(^^^)'Iq=( / {fiOyeU^,m~fi^,u)TK{{0,u),y)du-{y) 

v+.„ 



<2( J f{ey.,u;My)K{{0,u),y)du-{y) 
+ 2/ / p(u,u)K((ti,u),y)di/'(y) 



SO V is in fact a map from L'^{&) to L^(D'T). 

Then, following [1], we denote by Ly (D^) the closure of the set of gradients 
of all functions from L^((3). We then consider the orthogonal decomposi- 
tion of L^(9T) into the "potential" and the "solenoidal" subspaces: L^(Dl) = 
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I/y(0^) © (Ly (D^))-*-. To characterize the solenoidal subspace (Ly(Dl))-'-, we 
introduce the divergence operator in the following way. For / : i— M, we 
have div f : & ^ defined by 



(div/)(c.,w)= / K{iO,u),y)fiuj,u,y)du-iy) 

K{{0,u),y)f{ey.eU,Uy,{\ye\,u))du'^{y) (28) 



(note that for y G V^^^ we have {\y ■ e\,u) G ^LuiAy)- Now, we verify the 
integration by parts formula: for any / G L^((5), g G L^(D^) 

(^,V4^ = -(/div,>^. (29) 

Indeed, we have 

(^7,VA =( / K{{Q,u),y)g{uj,u,y){f{ey.,uMy)-f{^.u))du-{y) 



' KQ 

V, 



^((0, u), y)g{uj, u, y)f{uj, u) di^'^iy) 



+ ( j K{{Q,u),y)g{uo,u,y)f{dy.,uoMy)di^^{y)) . (30) 
v+,„ 

For the second term in the right-hand side of (|30|1 . we obtain 
^((0, m), y)g{uj, u, y)f{9y.eUJ,Uy) du'^{y)) 



d¥ da dfiQ^^{u,v)(^{0,u,a,v)g(uj,u,{a,v))f{9aio,v) 



1 

Z 

n A2 



— da dF dfi'^Q{v,u)^{a,v,0,u)g{eaUj,v,{\a\,u))f{uj,u) 



Z 



-oo f7 A2 
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f{uj,u)g{ey.eUjMy, i\y ■ e\,u))K{{0,u),y) di^'^iy) 



v: 

and so 



+ ( j f{iJ,u)g{9y.,u;,Uy,{\ye\,u))K{iO,u),y)du'^{y) 

and (!29|) is proved. 

Analogously to Lemma 4.2 of [4j, we can characterize the space (Ly (OT))-*- 
as follows: 

Lemma 3.4 g G if and only if div g{uj,u) = for Q-almost all 

{uj,u). 

Proof. This is a direct consequence of (!29l) . □ 

A function / G L^(DT) can be interpreted as a flow by putting formally 

f{uj,u,y) := -f{ey.eUJ,Uy, {\y-e\,u)) 

for ?/ G V~^, G ©. Then, it is straightforward to obtain that every / G 
Ly(D^) is curl-free, which means that for any loop yo,yi, . . . ,yn G du with 
yo = yn and K{yi, yi+i) > for i = 1, . . . , n - 1, we have 

n-l 

^f{ey^.eUj,Uyuyi+i - {yi ■ e)e) = o. (3i) 

Every curl-free function / can be integrated into a unique function : OJt t— *• 
M, which can be defined by 

n-l 

V?(u;, u,y) = ^f {Oy^.eUjMyh yt+i - (yi ■ e)e) (32) 

i=0 
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where yo, yi, ...,?/„ G So; is an arbitrary path such that yo = (0, u), yn = y, 
and K{yi,yi^i) > for 2 = 1, . . . , n — 1. Automatically, such a function ip 
satisfies the shift- covariance property: for any u G duo, x,y E du, 

ip{uj,u,x) = (p{uj,u,y) + ip[6y.eUJ,Uy,x - {y ■ e)e). (33) 

We denote by Ti.(JM) the set of all shift- covariant functions 9Jt M. Note 
that, by taking x = y = {0,u) in (133|) . we obtain 

(0,m))=0 for any ip G 7Y(0Jl). (34) 

Also, for any ip G H(DJl), we define gradip as the unique function / : — M, 
from the shifts of which (p is assembled (as in fl32p ). In view of fl34p . we can 
write 

{gradf){uj,u,y) = f{uj,u,y) for {uj,u,y) e%fe n{DJl). 

Let us define an operator C which transfers a function : 971 t-^ R to a 
function / : (5 i-^ R, / = C(p with 

{Cip){uj, u) = J K[{0, u), y) [(p{uj, u, y) - ip{u, u, (0, u))] dv'^{y). (35) 

Note that, by ([27]), we obtain C{V f) = Gf - f for any / G L'^{&). Then, 
using dMD and ([3lD, we write for ip G n{m) 

(divgradv3)(u;,u) = J K{{0,u),y)Lp{uj,u,y) du'^{y) 
v+ 

- J K{{0,u),y)Lp{ey.eUJ,l(y,{\y ■ e\,u)) du'^iy) 
K{{0,u),y)ipiiu,u,y) du^iy). 

Vj,„UVJ,„ 

So, for any (p G H(9Jl), we have divgradip = dp. This observation together 
with Lemma [3.41 immediately implies the following fact: 

Lemma 3.5 Suppose that ip G HiTl) is such that grady? G (Ly (D'T))-'-. 
Then (p is harmonic for the Knudsen random walk, i.e., {C(p){uj,u) = 
for Q-almost all {u!,u). 




24 



Now, we are able to construct the corrector. Consider the function 
(p{u!,u,y) = y ■ e and observe that ip G H{n). Let if = gradip. First, 
let us show that 

(^^>,, = li% (36) 
i.e., if (II2]) holds, then e L'^{m). Indeed, 



1 



+ 00 

2 



a da d¥ dfi^ ^{u,v)^{0,u, a,v) 



n A2 





^ J a'^ da J dF J dij!^ Q{v,u)^{a,v,0,u) 

-oo n A2 

/ {yeyK{{0,u),y)du-{y)), (37) 



V\^=l{ I {yerK{{0,u),y)du-{y)) 



and so 

Then, let g be the orthogonal projection of (— onto L|7(DT). Define 
ip e 7i(07l) to be the unique function such that g = gra.dip] in particular, 
ip[uj, u, (0, u)) = for M G duQ. Then, we have 

+ g = gTRd{{ye) + tlj{uj,u,y)) G (L^(O^))^, 

so Lemma [3.51 implies that for Q-a.a. {uj,u), ip is the corrector in the sense 
that 

E.((^i-^o)-e + ^(a;,n,ei)-^(^,M,eo) Uo = (0,^/)) =0 (38) 

(recall that, by ( !34l) . the term iP{ijj,u,^q) can be dropped from ( !38l) ). Now, 
denote 

J,- = E^((ei - ^o) ■ e + ilj{e^,.eu: Mio.il - {io ■ e)e)) \i^ = x). 
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By the translation invariance of P, fl38l) , and ([2]) , we can write 

+00 



A 



-co fl A 

1 

and this imphes that = for v'^ ® P-a.e. {uo^x). From (!33l) . we have 

ij{u,u,y) - ij{u,u,x) = ip{9^.eUJ,Ux,y - {x ■ e)e), 
which does not depend on u. We summarize this in 
Proposition 3.1 For ¥-almost all 00, it holds 

E.^ ((6 - ^0) ■ e + ^piuj, M, ^1) - ip{uj, u, ^0) \^o = x) =0, (39) 
for all u G dujQ and v'^ -almost all x G du. 

3.4 Sequence of reference points and properties of the 
corrector 

Let X be a random variable with uniform distribution in [—1/2, 1/2], inde- 
pendent of everything. Note that (x + ^, G Z) is then a stationary point 
process on the real line. For a fixed environment uo define the sequence of 
conditionally independent random variables ^ A, n G Z, with Cn uniformly 
distributed on du^ 



■'X+n- 



We denote by the expectation with respect to C and x (with fixed tu), 
and by E'' the expectation with respect to C conditioned on {x = 0}. Then, 
by (133|) we have the following decomposition 



n-l 



^(6'x^,Co, (^,Cn)) = ^^/'(6'^+iu;,0, (1,0+1)), (41) 



1=0 
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so that ipiO^uj, iP'i Cn)) IS a partial sum of a stationary ergodic sequence. 

By Condition L, there exists 71 > such that n'^^i^dun) > 71 P-a.s. Hence, 
since P is stationary and /€o,n G [0, 1], we obtain for / > 0, 



(E^/(e,^,Co)X = (E^/(^,Co)>, 



n A 

which imphes that 

if / G L\Q) then (E^/2( V, Co)>^ < 00. 

To proceed, we need to show that the random tube satisfies a uniform 
local D5bhn condition. Denote := K + K^^) + . . . + K^n) _ 

Lemma 3.6 Under Condition P, there exist N and 7 > such that for all 
x,y E TZui with \ (x — y) ■ e\ < 2 it holds that K^^\x,y) > 7, P-a.s. 

Proof of Lemma \ 3.6\ Indeed, with the notations used in Condition P and 
n = N -1. 



K^''\x,y) >K^'^^^\x,y) 

> [ K{x,zi)dv'^{zi) 



Bi 

X J K{zi,Z2) du'^{z2) . . . J K{Zn-l,Zn)K{Zn,y) dv"^ {Zn) 
B2 B„ 

□ 

Next, we state some integrability and centering properties, which will be 
needed later. 
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Lemma 3.7 We have 

{E^^'{u,u,{x,Co))\<'^, (43) 

(E«^(^^^,Co,(l,Ci))>, = (E^^(^,Co,(l,Ci))>, = 0. (44) 
Proof of Lemma \3. 7} We start proving that 

We know that grad-?/' G L^(DT), that is, ((grad?/')^) < oo. Analogously 

to the proof of the symmetry relation fl37|) we obtain (note also that, by fl29|) . 
(div^)^ = for all g e L\^)) 

((grad^)2>_ = ( / ^'iu,u,y)K{iO,u),y)du-iy) 
(E^(^2(a;,n,6) Uo = (0,m))> 



' KQ 



UJ.U 



1 

2* 



Then, using ( j3T|) we write 



n-l 

^^(w,«,^n) < n?/'^(u;,M,^i) + ra^?/'^(%^,.eu;,W^j,^j+i - {^j ■ e)e). 

Since Q is reversible for G, this implies that for any n 

( f ^'{cu,u,y)K^-\{0,u),y)du-{y)) < oo, (46) 



where K^"^^ (x, y) is the n-step transition density. 
Let us define 

= |x G (9cu : x-e G (n - l/2,n + 1/2]|, 

Now, we are going to use ( H6i) and Lemma [3.61 to prove f H5l) . Note that, by 
Condition L, there are positive constants Ci,C2 such that 

Ci < u'^iF^) < Ca, P-a.s. 
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Using (PT]) . we write on {x = 0} 

1 



7C1 

9aj 



+ / K^''^{il,Ci),yWiO,co,Ci,y-e)du-{y)). 



did 



Using the stationarity of C, under P, we obtain that 

(E^^2 (V, Co, (1, Ci)) = (E V Co, (1, Ci)) 
then, again by stationarity, 

/ ir(^)((l,Ci),Z/)^^(^i^,Ci,y-e)rfz/-(y) 



E< j K(^) ((0, Co), Co, y) du-{y))^. 

So, (gS]) follows from (H2l) and (Hel). 

Analogously, it is not difficult to prove that fl43l) holds. Indeed, similarly 



to (H2l) . we have 

1/2 

E^V^(t^,^i, (x,Co)) = / da—— — - dfi'^{v)ij'^{uj,u,{a,v)) 

-1/2 A 
1/2 

~ \ I / 

-1/2 A 

1 



7i 



i^,u,y)diy'^{y) 
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where we used in the last equahty. So, by Lemma IXBl 

111 J 

and (03]) follows from (|46D. 

Finally, let us prove (jS]). The first equality follows from the stationarity 
of P. Then, since gia-dip G Ly(D^), there is a sequence of functions e 
L^(6) such that V/„ grad?/^ in the sense of the L^(91)-convergence. Note 
that, in fact, when proving fHS]) . we proved that for any function g G 7i(07l) 
such that grad(7 G //^(S), we have for some C3 > 

(eV(V,Co,(1,Ci))>^ < C3((grad^)2)^^. 

Then, (jH]) follows from the above fact applied to g assembled from shifts of 
gradV' — V/„, since then we can then write 

(^(^^^,Co, (l,Ci))>, = £m [{Ue^^.u^Ql - (/„(e,.;,Co)>J = 

by the stationarity of P. □ 



4 Proofs of the main results 
4.1 Proof of Theorem 12.11 

In this section, we apply the machinery of Section [3] in order to prove the 
invariance principle for the (discrete time) motion of a single particle. 

Proof of Theorem \2.1[ Denote 

observe that by (15^ . is a martingale under the quenched law P^^. By shift- 
covariance fl33l) the increments of 0„ do not depend of x ^i-nd (. With the 
notation 

/i(a;,M) = E^[(0i-eo)' |eo = (0,M)], 
the bracket of the martingale 0„ is given by 

n-l 
1=0 
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By the Ergodic Theorem 

1 



-{QU-^a''^{hiuj,u)) (47) 

ft ^ 

a.s. as ^ oo. Clearly, o"^ G (0, oo). Moreover, for all e > 0, 



n-l 



-&^\> en'/' I e.] (48) 

for P-a.e. u and P(^-a.e. path. To show this, define for any a > and all 
n > 1 

hl:\iu) = E^((0„ - 0„_i)2l{|e„ - e„_i| > a} I a-i). 

Using the ergodicity of the process of the environment viewed from the par- 
ticle, we obtain 

1 " 

-E^^"^ - 0o)'l{|ei - Ool > a} I eo = (0,«))>^ 

i=l 

as n — > oo, for P-almost all u and P^-almost all trajectories of the walk. 
Note that, when a is replaced by en^/^, the left-hand side is, by Bienayme- 
Chebyshev inequality, an upper bound of the left-hand side of (HHI) multiplied 
by e^. Hence 0481) follows by taking a arbitrarily large. 

Combining fl47p and (1481) . we can apply the central limit theorem for 
martingales (cf. e.g. Theorem 7.7.4 of \TT\) to show that 



n 



-'/'ein-]^(TB{-) asn^oo, (49) 



where B{-) is the Brownian Motion. 

Then, the idea is the following: using (H^ and the Ergodic Theorem, we 
obtain that the corrector iI){uj,u,x) behaves sublinearly in x, which implies 
the convergence of n~'/''^[nt] ■ e. More precisely, we can write, with rrij : = 
[1/2 + ^j-e] and using 



©H ■ e + ^{dx^^ Co, ("^H, Cm[„t])) 

(50) 



^1/2 ^1/2 
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Let us prove that the second term in the right-hand side converges to in 
P(^-probabihty, for P-almost all oj and almost all ixX)- Suppose, for the 
sake of simplicity, that t = 1. Then, by the stationarity of the process 
((x + ra, Cn), n G Z) and f|T^ together with fHSl) . we have for alH > 

(E^E^ [V'(%-eo;,W6, {x + m-^^■ e, Cmj) | ^0 = (0, u)] )^ 

= (E^^2(a;,«,(x,Co))>^ 
< oo, 

and so, by the Ergodic Theorem 
1 " 

i=l 

as n OO, which implies that 

;^E^^^(6'x+Cn-e^>^^n, {rrinXm^)) (51) 

for P-almost all uo and almost all (x^O- Now, let us prove that the limit 
of the first term in the right-hand side of (150|) is the same as the limit of 

n~"^/^^[„t] ■ e; for this, we have to prove that 

> as n ^ oo, m Pa;-probabihty. (52) 



n 



1/2 



Using (jHj), flMj) . and the Ergodic Theorem, we obtain that for P-almost all oo 
m~^ip(6y^u}Xoi i^Xm)) for almost all {xX)y 1^1 ~^ oo. This means 
that, for any e > there exists H (depending on uj,(,x) such that 

\i/j{9^uj,Co,{m,(m))\ < H + e\m\. (53) 

Denote 

= ■ e + ip{Oy^uj, Co, ("^i, Cm,)) • 

From fl53l) we see that 
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^x^,Co, ("^j,Cm,))| < ^ + ^1 



<H+^- + e\ire\ 

<H+^- + e{\^j\ + \iP{e^ujXoAmjXm,))\): 
and so for e < 1/2 we obtain 

|^/'(^;,cu,Co,K-,Cm,))| < 2H + e + 2e\^j\. 
Using (jl9]) and ([5l]) in ([50]), we obtain 

|vl>.| 

max —7- a max |-B(s)|. 

i<n se[o,i] 



So, by the Portmanteau Theorem (cf. Theorem 2.1 (iii) of [2 

< P 



hm sup max \il)[e^uj, Co, {rrij, C^J) | > an^^^ 

n— >oo L i^:*^ 



max Ld s > — 

se[o,i] 2e 



which converges to for any a as e — > 0. This concludes the proof of Theo- 
rem EH □ 



4.2 On the finiteness of the second moment 

In this section, we prove the resuhs which concern the finiteness of (^)^- First, 
we present a (quite elementary) proof of Proposition 12. II in the case d> A. 
Proof of Proposition \2. 1\ (case d > 4). First of all, note that 

\{s G E,'^-^ : x + hs eRx A}\= 0{h-'-'^-^^) as /i ^ 00, 

uniformly in x G M x A. So, since C M x A, there is a constant Ci > 0, 
depending only on M = diam(A)/2 and the dimension, such that for P-almost 
all u 

P4\i^i -^o)-e\>h\^o = x]< C^h-^''-'\ (54) 

for all X G duj, h > 1. Inequality fl5^ immediately implies that b is uniformly 
bounded for > 4. □ 

Unfortunately, the above proof does not work in the case d = 3. To treat 
this case, we need some results concerning induced chords which in some 
sense generalize Theorems 2.7 and 2.8 of [B]. So, the rest of this section 
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is organized as follows. After introducing some notations, we prove Propo- 
sition I4.H which is a generalization of the result about the induced chord 
in a convex subdomain (Theorem 2.7 of [6]). This will allow us to prove 
Proposition 12.21 Then, using Theorem 2.8 of |6] (the result about induced 
chords in a general subdomain) we prove Proposition 14.21 — a property of 
random chords induced in a smaller random tube by a random chord in a 
bigger random tube. This last result will allow us to prove Proposition 12.11 
Let S" C A be an open convex set, and denote hj S = 'R x S the straight 
cylinder generated by S. Assuming that S C u, we let X be the event that 
the trajectory of the first jump (i.e., from to ^i) intersects S: 

I = {there exists t G [0, 1] such that + i^i ^ ^o)'t G S}. 

For any u G dS such that dS is differentiable in u, define n„ to be the normal 
vector with respect to dS at the point (0, u); clearly, we have n„ - e = (if dS 
is not differentiable in u, define n„ arbitrarily). Fix some family {Uy,v G dS) 
of unitary linear operators with the property UyC = hy for all v G dS. Now, 
on the event X we may define the conditional law of intersection of dS. 
Namely, for x, ?/ G du let 

t^^y = mf{t G [0, 1] : x + (y - x)t G dS}, (55) 

with the convention inf = oo. Then, we define the (projected) location of 
the crossing of dS by 



Uix + {y - X%^y), if e [0, 1], 

OO, otherwise, 



L{x,y) -- 

and the relative direction of the crossing by 
Y{x,y) 



0, otherwise. 



see Figure El Here, in the case when there is no intersection, for formal 
reasons we still assign values for L and Y; note, however, that in the case 
tx,y £ [0, 1] we have \-{x,y) G dS and Y{x,y) G §e- 

Before proving Proposition 12. 2[ we obtain a remarkable fact which is 
closely related to the invariance properties of random chords, cf. Theorems 2.7 
and 2.8 of [6]. We have that, conditioned on I, the annealed law of the pair 
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Figure 3: On the definition of L and Y: we liave L(a;, y) = Uz and ^{x, y) = 

of random variables (L(^0) ^i); Y(^0) ^i)) can be described as follows: L(^o! ^i) 
and y^io.ii) are independent, L(^O)^i) is uniform on dS and Y(,^0)^i) has 
the cosine distribution. More precisely, we formulate and prove the following 
result. 



Proposition 4.1 Letd>2. It holds that {P^[I]) = \dS\/Z. Moreover, for 
any measurable Bi C dS, B2 C Sg we have 

;E4l{L(eo,ei) G 5i, Y(eo,ei) G B^} I ^0 = (0,n))>^ = ^^7d j h-edh. 

(56) 



B2 



Proof. First, we prove (1561) . Define = {x G du : x ■ e G [—s, s]} for s > 0. 
By the translation invariance and ([2]), we have 

(E4l{L(eo,6) G 5i,Y(eo,6) G i?2} I eo = (0,n)))^ 
1 

~ Z 



dF J dfi^iu)K,l J du-iy)K{iO,u),y) 

n A dw 



xI{L((0,M),y)GSi,Y((0,M),2/)Gi?2} 



2sZ 



dFjdsj dfi^{u)Kll j du^{y)K{{s,u),y) 

Q -s A dui 
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2sZ 



xI{L((s,M),t/)G5i,Y((s,M),|/)G52} 
/ dv^{x) I du'^{y)I{l{x,y) G B,,Y{x,y) G B2}K{x,y). 



(57) 

Define the domain by 

V'^ = {x e u : X ■ e e [s, s]}, 

note that dV^ = U {{-s} x U ({s} x tu,). For x,y e dV^ let y) 
be defined as in ([3]), but with instead of u. Note that K{x, y) = K{x, y) 
when x,y E Ff . 

Next, we show that the random chord in uo with the first point in 
has roughly the same law as the random chord in V^: for any e > there 
exists So such that for all s > Sq (with some abuse of notation, we still denote 
by u'^l&D^) the {d — l)-dimensional Hausdorff measure of &D^) 



dp^{x) J di^'^iy)I{lix,y) G i?i,Y(x,y) G B2}K{x,y) 

du) 

dv^{x)dv-{y)l{L{x,y) G 5i,Y(x,y) G B2}k{x,y) 



< e 



(58) 



(in the second term, we suppose that L(x, ?/) = oo,Y(x, y) = when x G 
({— s} y. S)\J ({s} X S)). Indeed, we have 

z/"(F,") < z/"(9P^) < z/"(F,") + 2|A|, (59) 

and, by Condition L, there exists C4 > such that 

v^{F^) > C4S P-a.s. (60) 



Also, since C M x A, for any e > there exists C5 > such that for all 
X E duj 

J K{x,y)du'^{y) < e P-a.s. (61) 

{yeduj:\{x-yye\>C5} 
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Now, (135]) follows from (1^ - (1^ and a coupling argument: choose the first 
point uniformly on dV^; with big probability, it will fall on F^_c^ (and so 
it can be used as the first point of the random chord in du). Then, choose 
the second point according to the cosine law; by fl6T|) . with big probability 
it will belong to F^, and so the probabihty that the coupling is successful 
converges to 1 as s ^ oo. 

Then, recall Theorem 2.7 from [B]: in a finite domain, the induced random 
chord in a convex subdomain has the same uniform x cosine law. So, 

diy-{x)du^{y)I{L{x,y) G 5i,Y(x,y) G B2}K{x,y) 



I B I 

Pc^[^s]7^7d / h-edh, 



B2 

where X, is the event that the random chord of dV^ crosses the set [— s, s\ x 
dS. By the formula (47) of ^ (see also formula (17) in Theorem 2.8 there), 
we have 

2^^ (62) 

Since, by the Ergodic Theorem, |F^'^|/(2s) ^ Z as s — > oo, ( l62l) implies 
that P^[X,] \dS\/ Z a.s s oo. We obtain ([56]) using ([57]) and ([58]), and 
sending s to oo. 

Finally, the fact that (P<^[X]> = \dS\/ Z follows from ([56]) (take Bi = 
dS,B2 = Se). □ 

Now, using Proposition l4.11 it is straightforward to obtain Proposition l2.21 

Proof of Proposition \2.2i Suppose that u contains an infinite straight cylin- 
der S (more precisely, a strip, since we are considering the case d = 2) of 
height r > 0, P-a.s. Keep the notation tx,y from ( 155]) . and define also 

t'^^y = sup{t G [0, 1] : x + (y - x)t G OS}. 

On the event T, define the random points Tq, Ti G 95* by 

"^0 = ^0 + (^1 ~ ^o)tx,y, 
'^1=^0 + (^1 — ^o)t'x,y^ 
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so that (To, Ti) is the random chord of S induced by the first crossing. On X^, 
define Tq = Ti = 0. By Proposition l4.lt conditioned on X, the random chord 
(To, Ti) has the cosine law, i.e., the density of a direction is proportional to 
the cosine of the angle between this direction and the normal vector (which, 
in this case, is perpendicular to e). Let P[-] := ^{|P^[-I{X}])^ be the annealed 
probability conditioned on the intersection; since d = 2 and is a bounded 
interval, \dS\ = 2. With r] := (,^i — ^o)/||Ci ~ Co|| and n the inner normal 
vector to the cylinder at Tq, we have (see Figure HI) 



P[|(Ti-To)-e| >x] =P 



7] ■ n < 

Tv/2 



COS z dz 



X 



X 



2' 



and so the conditional density of the random variable | (Ti — Tq) - el is f{x) 
on M"*". Then, we have 



{\ = (E.(lte-eo)-eni{J} + I{n]|^o = (0,n))>^ 

> (E4|(6-^o)-epI{J}|^o = (0,n))>^ 

> (Ej(Ti-To)-en x(P.[J]) 
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Figure 5: Random chords induced in a random tube by a random chord 
in a random tube uj' (in this particular case, we have 6 = 4) 



Let us observe that if a stationary ergodic random tube is almost surely 
convex, then necessarily it has the form M x 5* for some convex (and nonran- 
dom) set S* C A. This shows that Proposition 14.11 is indeed a generalization 
of Theorem 2.7 of Now, our goal is to obtain an analogue of a more 
general Theorem 2.8 of [0]. For this we consider a pair of stationary ergodic 
random tubes {oj^uj') G VL^, let P be their joint law and P, P' be the corre- 
sponding marginals. Suppose also that u is contained in uj' P-a.s. We keep 
the notations such as Kx-, K{x, y), ■ ■ ■ for x,y & du' as well, when it can cre- 
ate no confusion; for the measures /i and v we usually indicate in the upper 
index whether they refer to uj or uj' . Denote also Z' = dF' K,oudfiQ {u). 

If (^0) ^i) is a chord in uj', we denote by (eS'\ • • • , the induced 

random chords in uj (see Figure E]). Here, i G {0, 1, 2, . . .} is a random vari- 
able which denotes the number of induced chords in uj, so that l = when 
the chord {^'o,^'i) has no intersection with uj. 

The generalization of Theorem 2.8 of [6] that we want to obtain is the 
following fact: 








which concludes the proof of Proposition 12. 2[ 



□ 
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Proposition 4.2 For any bounded function f : Wl t-^ W we have 
(e4/(^,W^o,6) Uo = (0,n))>^ 



f72 

X 



E..-' ( E / (\«.e^' - (ei'^ ■ e)e) I = (0, u)) . 

k=l 

(63) 

Proof. We keep the notations from the proof of Proposition 14.11 (with the ob- 
vious modifications in the case when u' is considered instead of a;). Without 
restriction of generality, we suppose that / is nonnegative. First, analogously 
to (1571) . we obtain that the right-hand side of (IMj) may be rewritten as 

2sZ J "-"(^'^'^ / du'^'ix) J du'^\y)K{x,y) 

^2 puj' dcu' 

^ S 

X E f{0.^^).e^Mx^'\ y^'^ - ix^'^ ■ e)e) =: Ti, (64) 

k=l 

where (x*-^-*, y^^-*), . . . , (x*^^*^^'^^-*, y*-'^"^'^-'-*) are the chords induced in cu by the 
chord (x, y) in u'. 

Let us denote Fj^,^, = {x E dcu : x-e E [h, h']} (so that F^ = F^^^J. Define 
t„(x, y) as the number of intersections of the chord (x, y) with F^_n s-n+i- To 
proceed, we need the following fact: let A be a subset of dV"^ and x E dV"^' . 
Then we have 

P(^/[random chord beginning at x crosses A\ < i^'^iA) sup K{x,y). 

Also, by decomposing A into pieces that may have at most one intersection 
with the chord starting from x and using the above inequality, we obtain 

Ei^i [number of intersections of the random chord from x with A] 

<u'^{A)snpk{x,y). (65) 

y€A 
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Using Condition L one obtains that z/'^(-F^„ s-n+i) bounded from above by 
a constant (see the argument before fl42p ). From we know that K{z, z') < 



Id 



\\z-z'\ 



7|T, so for any x G {s} x lu'^ it is straightforward to obtain that 



Ux,y)K{x,y)du-'{y) < ^^'^"(-^--"■--+1) < 9l, (66) 

n n 

Suppose without restriction of generahty that s is an integer number. Since 
/-{x, y) < 1 + X]n=i '^n{x, y), we obtain that 



-JdF{uj,J) j dv^\x) j dv-\y)k{x,y)i{x,y)<-^{l + Y,^ 



< (67) 

s 

and the same bound also holds if we change {s} x uj'^ to {— s} x uj'_^ in the 
second integral above. 

Note that, by the Ergodic Theorem, we have that 

^ ^ ^-Xavi-) ^ ^, ^^^^^^ p.^^^^ 

2s 2s 



Then, analogously to (l58l) . using fl67j) together with the fact that / is a 
bounded function, we obtain that for any e > there exists Sq such that for 
all s>so (recall ([MD) 

T2-T^< e, (68) 

where 

T2 ■■= -TT^wT / du-'ix)diy-'iy)k{x, y) 



X E /(^.«.e^,Wx(^'\2/^'^ - i^^'^ ■ e)e). (69) 

k=l 

Then, by Theorem 2.8 of [6], we have 

^2 = ^.(gp^) / du-ix)du-{y)k{x,y)f{e,.,u,Ux,y-ix-e)e). (70) 
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Again, it is straightforward to obtain that for any e > there exists Sq such 
that for all s > Sq 



1 



1 



du'^i^x) du'^{y)k{x,y)f{9^.eUJ,Ux,y- (x ■ e)e) 



2sZ 



dv-{x) dv-{y)K{x,y)f{e,.eUjMx,y-{x-e)e) < e. (71) 



lim / dv'^{x) I du'^{y)K{x,y)f{6^.eUJ,Ux,y - {x ■ e)e) 



By the Ergodic Theorem, we have that P-a.s. 
1 

s-^oo 2sZ 

fpi^ du) 
^ s 

= (E4/(cu,weo,ei)ieo = (o,M))>^, 

so, using (168|) . ( ITOi) . and ( ITTj) . we obtain, abbreviating for a moment 

fc=i 

that 

^2 A 

X E^,^,(2t I = (0,m)). 

(72) 

The other inequality is much easier to obtain. Fix an arbitrary c > 0, and 
consider 2ll{2t < c} instead of 2t. Since 21I{21 < c} is bounded, we now 
have no difficulties relating the integrals on F"^ x du' to the corresponding 
integrals on {&D^ )^. In this way we obtain that for any c, 

(E4/(^,W^o,ei) I ^0 = (0,n))>^ >^^L j dnu,u') j d^i-^{u)K,l^ 

n2 A 
xE^,,,(2lI{2l<c} |e; = (0,M)). 



We use now the monotone convergence theorem and f[72|) to conclude the 
proof of Proposition 14.21 □ 
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Using Proposition 14.21 we are now able to prove Proposition 12.11 for all 
d>3. 

Proof of Proposition lKli We apply Proposition 14. 21 with u' being the straight 
cylinder: cj' = M x A. For the random chord in a straight tube, using the fact 
that the cosine density vanishes in the directions orthogonal to the normal 
vector, we obtain that (for any starting point ^q) ipo := Ej^/((^[ — ^q) -e)^ < oo. 
Now consider the function fc{uj,u,y) = {y ■ eYl{{y ■ e)^ < c}. Since 

- ■ e)e) Q ■ e)^ 

fc=i 

we obtain that for any c 

(E4/,(a;,Weo,6) I ^0 = (0,n))>^ < ipo. 

Using the monotone convergence theorem, we conclude the proof of Propo- 
sition O □ 

Remarks, (i) Observe from the definitions of ipo above and ([1]) of A that 
ipo{M)=ipQ = M'^ipo{l). Then, we have shown the universal bound 

(b)^ < M'Cid), 

for aU^random tubes with vertical section included in the sphere A of ra- 
dius M, where C{d) = ipoi^) corresponds to the straight cylinder with spher- 
ical section of radius M = 1. 
(ii) For k > 1, denote by 

b^'\u,u) = EM^,-^o)-e\'\^o = iO,u)) 

the kth absolute moment of the projection of the first jump ]to the horizontal 
direction. Then, similarly to the proof of ]Propositions [2?T] and 12.2] one can 
obtain, for the rf-dimensional model, that (b^'^'') = +oo in the case when the 

random tube contains a straight cylinder, and that (6'-'^"'^)) < oo for any 
6>0. 
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4.3 Proof of Theorem [272] 

We start by obtaining a formula for the mean length of the first jump, at 
equilibrium. 

Lemma 4.1 We have 

i/2p/d+iw 1 r 
(E^llCi - U\\ = X i J l-ol dF. (73) 

Proof. Recall the notations F^, V^, K{x, y) from the proof of Proposi- 
tion I4.1[ The strategy of the proof will be similar to that of the proof of 
Proposition 14.21 Analogously to (1571) . we write 

U\ii-io\\\ = ^ JdrJ rf/i^(«)«:oi / di.^iy)K{{0,u),y)\mu)~y\\ 

n A duj 

1 



2sZ 



-s A 

X / du^iy)K{{s,u),y)\\{s,u)-y\ 



2sZ 

By Theorem 2.6 of [6j, we know that 



rfP / du-{x) / dv^{y)\\y-x\\K{x,y). (74) 



9a) 



^^/^Y{^)d 

du-{x)diJ^{y)K{x,y)\\x-y\\= ^,,\\ IP^I- (75) 



r(f + 11 



Denote by = {— s} x u^s and -D,. = {s} x u;^ the left and right vertical 
pieces of dV'^, so that dV^ = UDgU D,.. From ([71]) we obtain (recall also 
that K{x, y) = K{y, x) for all x,y e dV^) 

{E4^,-^,\\)^>^jdF j du-{x)du-{y)\\y-x\\K{x,y) 
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>^yrfFl J du'^ix)du-iy)\\y-x\\Kix,y) 

-2 j du-{y)\\y-x\\k{x,y)]. 

Observe that (recall ([I])) for all x, ?/ G R x A it holds that ||x — y|| < |(a; 
y) ■ e| + 2M. So, by (15^ . there exists Ci > such that for all s we have 



dv'^{y)\\y-x\\K{x,y) < Cilns + 2M, 

(_D<.UDr)xai5^ 

and, using fITS]) . we obtain 

(E.||6-^o||>^>^y rfP(^;^^^|Pr|-Ciln.-2M)rfP. (76) 



Since, by the Ergodic Theorem, 

> J \uJo \ dF a.s., as s — > oo. 



2s' ^ 



sending s to oo we obtain from ( 176|) that 

(E.llei-eo||>^>^^^^^xl|KMP. (77) 

Now, fix c > and write analogously to fl74|) 

(Eji^i-eop{liei-eoii <£}>^ 

= ^ /"rfP /" rf^^"(a;) /" rfz/-(y)||y-a;||I{||y-x|| <cK(a;,y). 



In this situation ||^i — ,^oP{||'Ci ~^o|| < c} is bounded. So, analogously to 
the proof of Proposition 14. II and again using fl75|) . by a coupling argument it 
is elementary to obtain that for any c 

l/2p/d±lW . p 

<E.||ei - eop{||6 - eoll < 5})^ < X - y |.;o| c^P. 
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Using the monotone convergence theorem and fl77p , we conclude the proof of 
Lemma 14.11 □ 

With Lemma 14.11 to hand, we are now ready to prove Theorem 12. 2[ 
Proof of Theorem \2.2 . Define n{t) = max{n : t„ < t}. We have 

r'/'Xt ■ e = r^'Hnit) ■ e + t-^/\Xt - ^(i)) ■ e. 

Let us prove first that the second term goes to 0. Indeed, by definition of 
the continuous-time process we have 

t-\{X, - ■ ef < ^{iUt)+i - Uit)) ■ e)'. (78) 

But then from the stationarity of ^ we obtain that 

n"^E^((^„+i - ^n) ■ e)^ -> 

as n — > oo for P-almost all u (this is analogous to the derivation of ( 15T|) in 
the proof of Theorem 12.11) . 

Now, the first term in the right-hand side of (ITHI) equals 



n(t)y/2 1 ^ 



t J n{t) 

For the second term in the above product, apply Theorem 12. 1[ As for the 
first term, since 

<t) ^n(fy^ n{t) 



Tn(t)+1 t Tn{t) 

by the Ergodic Theorem and Lemma 14.11 we have, almost surely, 
lim — ^ = lim — 

t— >oo t n^oo Tn 



((E^ll^i-^olDJ 

\uJo\d¥ 



>2 



7ri/2r(^^)d 

This concludes the proof of Theorem 12. 2[ □ 
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Appendix 



In this section we discuss the case when the map a ^— cJq, is not necessarily 
continuous, which corresponds to the case when the random tube can have 
vertical walls. The proofs contained here are given in a rather sketchy way, 
without going into much detail. 



to be the section of the boundary by the hyperplane where the first coordinate 
is equal to a. Then, let 



Besides Condition R, we have to assume something more. Namely, we 
assume that for P-almost all a;, v'^-alm.osi all (a, u) G duo are such that either 
\wa\ > (so that a E 5*^"^ for some n), or {a,u) E IZ^ (recall the definition 
of TZ^ from Section [2]). 

Also, we modify the definition of the measure /x^ in the following way: it 
is defined as in Section [2] when = 0, and we put /i^ = when \zUa\ > 0. 

Observe that, for any n > 0, iS^"-* is a stationary point process. Note that, 
in contrast, the set U„>oiS^"'^ may not be locally finite, which is the reason 
why we need to introduce a sequence Let P^"") be the Palm version of P 
with respect to iS'-"'\ i.e., intuitively it is P "conditioned on having a point 
of at the origin". Observe that P*^"^ is singular with respect to P, since 



Define 



Wa = 



{u E A : (a, u) E duj} 



5^°) = {a G M : > 1}, 



and, for n > 1 





)}) 
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for n > 1. Then, we define the measure Q which is the reversible measure 
for the environment seen from the particle: 



u},u) = — 



uar'- '[u 



n=0 



(79) 



where Z is the normalizing constant; as we will see below, Z still can be 
defined directly through the limit 



Z = lim 

a^oo 



'{{x e duj : X ■ e e [0,a]}) 



The scalar product is now defined by 



1 



/ duf{u,u)g{u,u) 



n A 

oo 

71=0 



Now, we need a slightly different definition for the transition density: define 
K{x,y) by the formula ([3]) but without the indicator functions that |n^(a;) • 
e| 7^ 1 and \n^{y) ■ e\ 1. Also, the transition operator G can be written in 
the following way: 

Gf{u;,u) = E^(/(%.e^,W^i) I ^0 = (0,n)) 

K{{0,u),x)f{e^.,uj,Ux)diy-{x) 



dui 



da j d^C{v)Kllf{e^uj,v)K{{f},u),{a,v)) 

-oo A 

oo „ 

+ E E dvf{e^uj,v)K{{0,u),{a,v)). 



n=0 Q:e5(")^_ 

Now, we have to prove the reversibility of G with respect to Q. The direct 
method adopted in the proof of Lemma 13.11 now seems to be to cumbersome 
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to apply, so we use another approach by taking advantage of stationarity. 
For two bounded functions /, g, consider the quantity 

21(a) = ^ / du^{x)du^{y)K{x,y)f{e^.,ujMx)g{ey.,ujMy)- 
jf,a J 

{ze9w:x-ee[0,a]}2 

Using ( 16T|) . it is elementary to obtain that (assuming for now that the limit 
exists P-a.s.) 

lim 2t(a) = lim ^ / rfz/"(x)/(^^.eCJ, Wx) 

a^co a^co Zd J 

{a;£9<^:a;-eS[0,a]} 

X j du'^iy) K{x,y)g{9y.eUJ,Uy) 

du) 

= lim ^ / du'^{x)f{9:r;.eUJ,Ux)Gg{9^.eUJ,Ux). 

a^co jf,a J 

{x£Oaj:a;-eS[0,a]} 

Then, we write using the Ergodic Theorem 

a 

lim - I da I djj!^{u) K~\f{6aUJ,u)Gg{6aUj,u) 

a— ►oo a J J 

A 

= jd^J dfx'^iu)K^lf{LU,u)Gg{LU,u) p-a.s. 

Q A 

Again, by the Ergodic Theorem, we have 

n l^^'") n [0,a]| 

lim J = Qm P-a.s., 



so that we can write 



lim- / du f{6aUJ,u)Gg{6ato,u) 

a^oo a J 
ae5{'")n[0,a] „^ 

= qm^ rfP^""^ y duf{uj,u)Gg(u:,u) P-a.s. 
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Thus, we have for P-almost all environments 



lim 2t(a) = lim — 

a— >oo a— »oo Zd 



da / dn'^{u) K^ J{9aU},u)Gg{9aUJ,u) 



A 

oo 



m=0 ae5(™)n[0,a] 



du f{OaUJ, u)Gg{9aLJ, u) 



1 



(iP / d/x^(n)/s:o,i/(a;,n)G'^/(^,n) 



n A 

oo 



m=0 n „ 



By symmetry, in the same way one proves that lima^oo 21(a) = ((7, Gf)^, and 
so G is still reversible with respect to Q. 

Now, the crucial observation is that the formula fH2l) is still valid even in 
the case when Q is defined by ( 17^ . since we still have, for any / > 



so one can see that all our discussion goes through in this general C3iSG clS well. 
However, we decided to write the proofs for the case of random tube with- 
out vertical walls to avoid complicating the calculations which are already 
quite involved. Here is the (incomplete) list of places that would require 
modifications (and strongly complicate the exposition): 

• the display after fl30|) (part of the proof of (l29l) ): 

• the proof of 

• the proof of Proposition 13.1} 

• the proof of ( H3i) : 

• calculations fl57|) and fITil) . 
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